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Recent results of rare B decay analyses based on 3f .9 million BB collected with the Belle 
detector at the KEKB asymmetric e + e~ collider are presented. We have made the first 
observation of charmless baryonic decay B — > ppK^ , the three-body B° 
B° — > K° K + K~ . The measured branching fractions are B(B 

W- 6 , B(B°^K°n + n-) = (53.2 ± 11.3 ± 9.7) x 10~ 6 , and B(B°—>K°K + K~) = (34.8 ± 
6.7 ± 6.5) x 10~ 6 . We also see strong evidence of B ± — > r/K ± and B ± — > rjn ± , and observed 



ppK+) = (4.3+ i ± 0.5) x 



the decay B -> wif* with B(B ± -> wA"*) = (9.9l2;I ± 1.0) x 10" b . Preliminary results 



of improved measurements of the branching fractions for the decays B - 
reported. No evidence for direct CP violation is found in the decays B 
K ^. ttV, and uo^. 



Kn and -kit are 



1 Introduction 

One of the most important goals of experiments at B-factories is to precisely measure the sides 
and angles of the unitarity triangle in the Cabbibo-Kobayashi-Maskawa matrix! (CKM) and 
check its consistency. Any inconsistency is a clear signal of new physics beyond the Standard 
Model (SM). Rare B decays play an important role in working towards this goal. 

Preliminary results of rare B decay analyses basecLon a 31.9 million BB sample are presented 
here. The data were collected with the Belle detectorlil at the KEKB asymmetric e + e~ collider. □ 

Belle is a general-purpose detector with a 1.5 T superconducting solenoid magnet. Charged 
particle tracking, covering 92% of the total center-of-mass (CM) solid angle, is provided by the 
Silicon Vertex Detector (SVD) consisting of three concentric layers of double-sided silicon strip 
detectors and a 50-layer Central Drift Chamber (CDC). Charged hadrons are distinguished by 
combining the responses from an array of Silica Aerogel Cerenkov Counters (ACC), a Time of 
Flight Counter system (TOF), and dE/dx measurements in the CDC. The combined response 
provides K/n separation of at least 2.5a for laboratory momenta up to 3.5 GeV/c. Photons and 



Table 1: The signal yield, reconstruction efficiency (e), statistical significance (£), branching fractions (B), and 
the 90% confidence level upper limits (UL) for two-body B decay modes. All the results are preliminary. 
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electrons are detected in an array of 8736 CsI(Tl) crystals (ECL) located inside the magnetic field 
and covering the entire solid angle of the charged particle tracking system. The 1.5 T magnetic 
field is returned via an iron yoke, instrumented to detect muons and Kl mesons (KLM). The 
KLM consists of alternating layers of resistive plate chambers and 4.7 cm thick steel plates. 

In all the decay modes presented here, the continuum process (e + e~ — > qq) is the dominant 
background. Since BB events are spherical while the continuum events are jet-like, we apply 
cuts on various event shape variables (such as sphericity, thrust angle, Fox- Wolfram moments, 
and the production angle of B) to suppress the background. 

B candidates are identified using two kinematic variables: beam constrained mass: M^ c = 
V E beam ~ pI> and tne energy difference: AE = E B - E b 

eam- Here E bearn is the beam energy, 
Pb and Eb are the momentum and energy of a reconstructed B candidate, respectively, where 
all variables are defined in the T(4S') rest frame. 

K/ir separation is performed by applying a cut on the likelihood ratio, Lk/(L 1t + Lk), where 
Lk (L n ) is a kaon (pion) likelihood computed from information from the particle identification 
devices: specific ionization loss in the central drift chamber, photo-electron yield in the aerogel 
Cherenkov counters, and time-of- flight. § 

2 Two-Body Charmless Hadronic B Decays 

These processes are manifestations of penguin or suppressed three amplitudes proportional to 
small couplings in hadronic flavor mixing (CKM matrix). Because of the absence of CKM favored 
b — > c amplitudes, these decays are particularly sensitive to potentially new contributions from 
interference effects and virtual particles in loops. 

2.1 B^7nr,Kir,KK 

The charmless hadronic B decays B — > -kit, Kit, and KK provide a rich sample to test the stan- 
dard model and to probe new physics. Of particular interest are indirect and direct CP violation 
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Figure 1: AE distributions of B — > K-k and rnv decays: (from left to right) K + ir , K + -k°, Kg-ir + , Kgir , tt + -k , 

tt+tt , and ttV. 



in the irir and Kir modes, which are related to the angles 02 an d 4>3 of the unitarity triangle, 
respectively. □ Measurements of branching fractions of these decay modes are an important step 
towards these CP violation studies. 

In this analysis, B meson candidates are reconstructed in ten combinations: h + h'~, /i + 7r°, 
K° s h+, K° s ir°, K° S K° S , and ^°vr°, where the symbols h and U refer to ir or K. Candidate it 
mesons are formed from pairs of photons and candidate Kg mesons are reconstructed from pairs 
of oppositely charged tracks with a displaced vertex from the interaction point. Signal yields 
are extracted by fits to the AE distributions taking into account feed-across from other mis- 
identified B — > hh' decays and backgrounds from multi-body and radiative charmless B decays. 
The fit results are given in Table |l] and the AE distributions are shown in Figure ||. 



2.2 B 



± 



Previous measurementsM yielded large rates for B — > r]'K and B — > r]K* , motivating a number 
of new theoretical ideas. Measurement of related decays B — > r\K can help clarify these. Besides, 
it has been suggested that the decays of B + — > r]ir + is a good candidate for observing direct CP 
violation. 

In this analysis, we reconstruct r/ mesons using the r/ — > 77 and r\ — > 7r + 7r~7r decay channels. 
Candidate r/ mesons are required to have invariant masses within ±2. 5a of the 77 peak, where a 
is 10.6 MeV/c 2 and 3.4 MeV/c 2 for the 77 and 7r + 7r~7r° modes, respectively. For the 7r + 7r~7r° 
mode, the 7r + 7r _ pair is constrained to a vertex. Both photons from the rj — > 77 mode are 
required to have E~ > 100 MeV and we remove rj candidates if either of the daughter photons 
can be combined with any other photon with Ej > lOOJVIeV to form a tt° candidate. The rj 
candidates are further constrained to the known 77 mass. 

Signal yields are obtained from extended unbinned maximum likelihood (ML) fits on the 
variables Mb c and AE. The signal probability density functions (PDF) are Gaussian for Mb c 
and an empirically determined parameterization □ for AE. The background PDF are taken to be 
an empirical function0 for Afb c and a first-order polynomial for AE. The statistical significance 
is defined as y/— 21n(£(0)/£ max ) where £ max is the likelihood at the nominal signal yield and 
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Figure 2: Mb c and AE distributions for (left) B — > T)K and (right) B — > ■r]-K ± . Histograms represent data, 
with the rj^+^-^o subset shaded, the solid curves represent the fit functions. 
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Figure 3: The Mj, c and AE projections with the normalized signal and background PDFs for (left) B^ 1 — > uiK 

and (right) B —* urn . 

C(0) is the likelihood with the signal yield fixed to zero. The results of the fits for yields are 
given in Table |l[ Figure |2| shows distributions of Mb c and AE with the projections of the fit 
function. We see strong evidence for the decays — > r\K^ and — > rjir^ with statistical 
significance of 4.9<7 and 4. 3a, respectively. Assuming the number of B°B° and B + B~ pairs 
to be equal, we find their branching fractions B{B^ — > rjK^) = (5.3l^"g ± 0.6) x 10 -6 and 
B(B ± — > 777r ± ) = (5.4+^ ± 0.6) x 10 -6 . 

2.3 B ± -> uh^ 

In this analysis, candidate to mesons are reconstructed from 7r + 7r~7r combinations where the CM 
momentum of the ir° is required to be greater than 350 MeV/c to reduce the large combinatorial 
background from low energy photons. The invariant mass of the 7r + 7r~7r° combination is required 
to be within ±30 MeV/c 2 of the nominal uj mass@ (the natural width of the uj is 8.9 MeV/c 2 ). 

Signal yields are extracted using unbinned ML fit simultaneously for M^ c and AE. The 
background functions include a combinatorial component and a component from other charm- 
less B decays. Due to the possible mis-identification between K~ and 7r~, we also include a 
component for feed-down from ujtt^ to loK~ fitting and vice versa. If a pion is mis-identified 
to a kaon, the wrong mass assignment shifts the ojit~ signal 44 MeV away from zero, and the 
feed-down component from luit~ can be distinguished in the AE distribution. The results of 
the fit are summarized in the Table |l[ The projections of the M& c and AE with the normalized 
signal and background PDFs are shown in Figure ||. We have observed the decay B^ — ► ojK^ 
with 6.4u significance. Our results on the B~ — > luK~ and u>tt~ branching fraction measure- 
ments disagree with the earlier ones from CLEO0 and BABArI But the sum of the branching 
fractions of B~ — > {ojK~ + uon~) is consistent with the CLEO's result. 

3 Three-Body Charmless Hadronic B Decays 

Belle has recently published results on tree body charmless hadronic decays B^ — > ir ± /t + /t _ El. 
It is of interest to look for similar phenomena in the neutral channel, to further investigate the 
b — ► s penguin transitions which mediate these decays. In addition, these modes may in future 
be used to measure CP violation. 

3.1 B° K h+h~ 

In this analysis, we reconstruct the decays B° — ► K°h + h~ without any assumption on the 
intermediate hadronic resonance. Only K° — > Kg — > TT + ir~ is considered here. As in other 
rare B decay modes, continuum events are the dominant background source and are suppressed 
using various event shape and kinematic variables. In the B° — > K°tt + tt~ mode, we also have 
large backgrounds from B®—>D~i: + , D~^>Kstt~ and B®—>J/i[)Ks, J where muons are 
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Figure 4: The AE (left) and Mb c (right) distributions for three-body final states: Kstv + tv (top), KsK + K 

(middle) and Ka^if* (bottom). 

misidentified as pions. These backgrounds are suppressed by requiring \M(Ks^) — M D - \ > 
0.100 GeV/c 2 , \M(h+h~) - Mj/^\ > 0.070 GeV/c 2 , and \M(h + hr) - M^ [2S) \ > 0.050 GeV/c 2 , 
where h + and h~ are pion candidates. Signal yields are obtained from fits to the AE distribu- 
tions. We find 60.3 ±11.0 B° -» K°ir + ir- events and 57.9 ±10.0 B° -> K°K + K~ events, which 
corresponds to preliminary branching fractions: 

BOB -^ ^ - ) = (53.2 ± 11.3 ± 9.7) x 10" 6 , (6.6a) 
B{B°^K°K + K-) = (34.8 ± 6.7 ± 6.5) x 10" 6 (7.4a) 

Figure || shows the AE and Mt, c distributions for these decays after the background suppression. 
We also search for B° — > K^K^-k^ but do not observe significant signals. The 90% C.L. braching 
fraction upper limit is calculated to be B(B° -» iT ^^) < 9.3 x 10~ 6 . 

Further studies of intermediate resonant states of these decays are made with a Dalitz plot 
style analysis. Figure |5| shows the 7r + 7r - and K^tt^ invariant mass distributions for selected 
B° — > Kg7T + 7T^ candidates in the B signal region, and the K + K~ and K^K^ invariant mass 
for B° -» K° S K + K-. Clear contributions from K*(892) ± vr =F and 0(1020)^° are seen. We 
also find broad resonances in Kgtr^ and K + K~ mass around 1.4 GeV/c 2 and 1.5 GeV/c 2 , 
respectively. Assuming a set of two-body final states, we perform a simultaneous likelihood fit 
to the AE distributions for various resonance regions and determine the exclusive branching 
fractions The uncertainty due to possible interference between different intermediate states is 
included as a model-dependent error. The results are summarized in Table |2[ 
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Figure 5: (a) 7r + 7r and (b) 7r + Kg invariant mass spectra for selected B° — > K%tv + ti , (c) K + K and (d) K + K% 
for B° — > KgK + K~ , signal events (open histograms) and for background events in the AE sidebands (hatched). 



Table 2: Preliminary results of the simultaneous fits to the B — > K tv n and K K K final states, respectively. 
Branching fractions of the corresponding decay modes from B + — » K + h + h~ are also listed for comparison. 13 
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4 Charmless Baryonic 5 Decays 

In contrast to charm meson decay, final states with baryons are allowed in B meson decay. To 
date, a few low multiplicity B decay modes with baryons in the final state from b — > c transitions 
have been observed. c3 Rare B decays due to charmless b — > s and b — > u transitions should also 
lead to final states with baryons. A number of searches for such modes have been carried out by 
CLEO,0 ARGUS, 111 and LEP0 but only upper limits were obtained. Stringent upperlimits 
for two-body modes such as B — > pp, Ap, and AA have recently been reported by Belle. E-3 

4.1 B ± ^ppK ± ^ 

We have searched for the decay modes B^ — ► ppK^ and B° — ► ppKg. These modes are expected 
to proceed mainly via b — > s penguin diagrams. We also search for B + — > ppir + which is expected 
to occur primarily via a b — > u tree process. Once they are established, these baryonic modes 
may be used to either constrain or observe direct CP violation in B decay. 

To reconstruct signal candidates in the B + — > ppK + mode, we form combinations of a 
kaon, proton and anti-proton that are inconsistent with the following b — > ccs transitions: 
B + -> J/ipK + , J/ip -» pp; 5+ -> r/ c A"+, r/ c -> pp; -> ^' -> PP and 5+ -> x c [o,i]^ + , 

Xc[o,i] — * PP- This set of requirements is referred to as the charm veto. Similar charm vetoes are 
applied in the analysis of the other decay modes. In the case of B° — > ppK$, events with pK$ 
or pKg masses consistent with the A c are rejected. 

In Figure §, we show the AE distribution (with 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 ) and 
Mb c distribution (with \AE\ < 50 MeV) for the signal candidates. We fit the AE distribution 
with a double Gaussian for signal and a linear background function with slope determined from 
the Mb c sideband. The fit to the AE distribution gives a yield of 42.8^° 6 8 with a significance 
of 5.6a. This is the first observation of a b — > s transition with baryons in the final state. 

We also examine the M p p mass distributions for events in the AE, M^ c signal region. The 
signal yield as a function of pp mass is shown in Figure ||(c). These yields were determined by fits 
to the AE distribution in bins of pp invariant mass. The distribution from a three-body phase 
space MC normalized to the area of the signal is superimposed. It is clear that the observed 
mass distribution is not consistent with three-body phase space but instead is peaked at low pp 
mass. This feature is suggestive of quasi two-body decay. It is also possible that the decay is 
a gemiiri e t hree-body process and that this feature of the Mpp spectrum is a baryon form factor 
effect. HI 

To avoid model dependence in the determination of the branching fraction for ppK + , we fit 
the AE signal yield in bins of M(_pp and correct for the detection efficiency in each bin using a 




Figure 6: (a) AE and (b) Mb c distributions for B + — > ppK + candidates, (c) The fitted yield divided by the bin 
size for B + — > ppK + as a function of M v p. The distribution from non-resonant B + — » ppK + MC simulation is 
superimposed. The inset shows the M p p distribution for the J/ipK + signal region by removing the charm veto. 

three-body phase space B + — > ppK + MC model. We then sum the partial branching fractions 
in each bin to obtain 

B(B + -^ppK + ) = (4.3+^(stat) ±0.5(syst)) x 1(T 6 . 

To verify the analysis procedure and branching fraction determination, we remove the J ftp 
veto and examine the decay chain B + — > J ' /tpK + , J/tp — > pp. The pp invariant mass spectrum 
for J/ipK + signal candidates is shown as an inset in Figure |6](c). The obtained branching 
fraction is in good agreement with the PDG world average. 

For B° — » ppKg and B + — * ppir + modes, after the application of the charm and A c ve- 
toes, no significant signals are observed. A fit to the AE distributions to ppK® gives 6. 4^3'^ 
events and toot 111 gives 16.2^'q events with significance of 2.1a. Applying the Feldman-Cousins 
procedure, c3, we obtain an upper limit at 90% C.L. of B(B° — > ppK°) < 7.2 x 10 -6 and 
B{B+ -> ppTr+) < 3.7 x 10" 6 . 

5 Search for Direct CP Violation 

The most straightforward indication for CP violation in the B meson system would be a time- 
independent rate asymmetry between CP conjugate decays into flavor specific or self-tagging 
final states. Direct CP violation (DCPV) of this type will occur in a decay containing at 
least two amplitudes that have different CP conserving and CP violating phases. The charge 
asymmetry will be most evident in decays where the two amplitudes are of comparable strength. 
The partial rate asymmetry can be written as 

__ N(B -> /) - N(B -» /) _ 2|^4i||^ 2 | sin5sin0 

CP ~ N(B -> /) + N(B -»■ /) ~ |Ai| 2 + |^ 2 | 2 + 2|A 1 ||A 2 | cos(5cos^' 

where 5 and <j> are the CP conserving and CP violating relative phases between amplitudes A\ 
and A2', B represents either a B^ or B + meson, / represents a self tagging final state, and B 
and / are the conjugate states. In the Standard Model, DCPV occurs in charmless hadronic 
B decay modes that involve both penguin (P) amplitudes and weak b — * u tree (T) amplitudes 
containing the CP violating weak phase <ps = arg (V* b ) (in a standard conventions). 

We search for direct CP violation in the B — > K + ir~ , K + ir°, Kgir + , and uK + decay modes 
by measuring the difference between the yields of B and B decays into the self tagging final 
states. The asymmetries obtained are summarized in Table ||. All the results are consistent 
with null asymmetry, hence we set 90% C.L. limits. Note that the systematic bias from charge 
asymmetry in the detectors is less than 1%, much smaller than the statistical errors at present. 



Table 3: Results of searches for direct CP violation. 



mode 


N(B) N(B) 


Ac P 90% C.L. 


K+7T- 

K+n° 
K° s Tr+ 
n+ir° 


103 ±12 115 ±14 
28 ± 8 30 ± 8 
49 ± 8 18 ± 6 
24 ± 8 13 ± 7 


-0.06 ±0.08 ±0.01 -0.20:0.09 
-0.04 ±0.19 ±0.03 -0.39:0.30 
0.46 ±0.15 ± 0.02 0.18:0.73 
0.31 ±0.31 ± 0.05 -0.25:0.89 


ujK+ 


7.4 ±3.5 11.6 ±3.7 


-0.22 ±0.27 ±0.04 -0.70:0.26 



6 Summary 

We have made the first observation of charmless baryonic decay — > ppK , the three-body 
B° — > K°7r + ir~ and B° — > K°K + K~ , and see strong evidence of -B 1 * 1 — > ry^T^ and £? — ► 777r =l= . 
We also observed the decay B^ — > ujK^ and measured the branching fractions for the decays 
B — > if7r and irir. We find no evidence for direct CP violation in the observed decays. By 
summer 2002 it is anticipated that Belle will have collected 90 fb _1 of data providing a rich 
sample to continue the search for rare B decays and measure CP violating effects in a variety 
of B decay modes. 
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